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Abstract. We derive the renormalization group equations for dimensionless couplings and soft 
supersymmetry breaking (SSB) parameters given arbitrary high scale inputs, incorporating 1-loop 
thresholds. A program (to be incorporated into ISAJET) has been developed, which includes 
flavour violating couplings of superparticles based on an arbitrary framework for the high scale 
flavour structure of SSB parameters. The weak scale flavour structure derived in this way can 
be applied to the study of SUSY flavour changing decays. As an illustration, we recalculate the 
branching ratio of the flavour-violating decay of the top squark. 

PACS. ll.30.Pb Supersymmetry - 11. 10. Hi Renormalization group evolution of parameters - 
14.80.Ly Supersymmetric partners of known particles 



1 Introduction 

In most supersymmetric models, the soft SUSY break- 
ing parameters are determined by high scale physics. 
Using an appropriate high scale ansatz, weak scale 
SUSY couplings relevant for phenomenology can be 
calculated using the renormalization group equations 
(RGEs). This procedure also yields predictions for flav- 
our violation at the weak scale, assuming of course we 
have a theory (or at least an ansatz) of flavour at a 
high energy scale. 

To accurately obtain the RG prediction at the two 
loop level, effects of thresholds at one loop must be 
incorporated. These include the decoupling of heavy 
particles and the concomitant splitting between cou- 
plings that are equal in the SUSY limit (the gauge and 
gaugino couplings, for instance). These new couplings 
will be identified with a tilde and labelled with the 
scalar which appears in the operator, e.g. the Buuji 
coupling is denoted by {g'u^)- 

The fact that gaugino couplings can become dis- 
tinct from their corresponding SM gauge counterparts 
means that they can even develop flavour off-diagonal 
terms. These couplings therefore play an important 
role in flavour physics. It is also important to remem- 
ber that threshold effects must be included in the run- 
ning of the soft masses and trilinear SSB parameters. 

The RGEs are constructed to describe a collec- 
tion of effective theories (with varying particle content) 
valid at different scales. We approximate the chang- 
ing particle content via step functions which remove 
the influence of the heavy particles at the appropriate 
scales lUE]- Each threshold is denoted hy a 6 and la- 
belled with the relevant particle. Above all thresholds, 
in the unbroken regime, all ^-functions are unity, and 
the running is as given by the MSSM |3j. Below all 
sparticle and additional Higgs thresholds, these step 
functions all vanish and the running is purely SM [4]. 



The importance of considering these effects can 
be shown by re-examining previous estimates for the 
flavour changing decay ti — > cZi, which requires a so- 
lution of the RGEs in order to estimate t — c mixing. 
This has previously been estimated [5] by integrating 
the RGEs in a single step. We will show that the results 
are significantly different when the RGEs are solved 
numerically, and the various effects discussed above 
are taken into account. 



2 Deriving and solving the threshold RGEs 

The general procedure for inserting thresholds begins 
with RGEs for a general theory [6] which do not de- 
pend on the particle structure of the theory. The RGEs 
in [6| are written for dimensionless parameters. In [7] 
they are rewritten for complex fermions and expanded 
to include RGEs for dimensionful parameters which 
are necessary for supersymmetric theories. We write 
the interaction of 2-componcnt Majorana fermions (tpp) 
with real scalars {4>a) via Yukawa couplings (Ip^) as 



Hi 



(1) 



using C as the metric which makes the fermion bilinear 
Lorentz invariant. With this interaction, the correct 
form of the general Yukawa RGE is 



dt 



1 — loop 

i [Y^(F)Y'^ Y"Y2(i^)] + 2Y''Y^"Y'' 
_l_Y''Tr 1^ — (y'^^Y'^ + Y^^" Y** 
-3g2{C2(F),Yn, (2) 
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in which Y2(F) = Y^''Y^ and C2(F) = f^t^. An ex- 
ample, rewritten in the more phenomenologist friendly 
4-component language with complex scalars and both 
Majorana and Dirac fermions, is shown in AppendixlAl 
It should be stressed that though the equations look 
much lengthier when written in 4-component form, 
much work needed to derive the RGEs has already 
been done, so their larger size is also indicative of ease 
of use. 

To facilitate decoupling, interactions of the fields 
need to be written in the (approximate) mass basis. In 
the Higgs sector the real neutral fields h, G", H and 
A are combined into two complex Higgs fields, h and 
Ti. (using a different font to differentiate the complex 
fields from the real fields) such that 



n 



h+iG° 
V2 
-H+iA 
V2 ■ 



(3) 

(4) 



h and H along with the complex charged fields 
and are obtained from the usual Higgs doublets 
(in the notation of 8J) through a field rotation: 



G+ 
h 

H+ 
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(5) 



(6) 



where s = sin /? and c = cos (3. We consider the ef- 
fects of Higgs boson thresholds in the approximation 
mA » Mz, in which case mn ''^ ruA ^ ^h+i with 
rufi ~ M^Q. The doublet in ([6]) then decouples at a 
scale Q mn while that in ([5|) remains in the low 
energy theory. 

For the neutralino and chargino thresholds, we ig- 
nore gaugino-Higgsino mixing (an excellent approx- 
imation if either of the conditions >> Mw or 
|Afi,2| >> are satisfiec0) so that the Higgsino 

mass eigenstates 



V2 



(7) 



are decoupled at the scale while the bino and wino 
states are decoupled Q — Mi and Q — M2 respec- 
tively. 

Below Q ^ rriH the theory contains just the SM 
Higgs scalar (although it may still contain Higgsinos), 
and the MSSM Yukawas (/) are replaced by their SM 
equivalents (A) using the conditions: A„ = sin/3/„ and 
•^{d,e} — cos/3/{rf.e}- Notc that if < ruH, Higgsinos 
would still couple fermions to sfermions, an example 
of which is the term Wi^ouun with coupling (/"^). 



If instead niA ~ Mz, all Higgs bosons have a roughly 
common mass, and threshold effects in this sector are unim- 
portant. 

^ If neither of these conditions are met, the gauginos and 
Higgsinos will all have similar mass so they will all decouple 
at about the same scale. It will therefore still be approxi- 
mately valid to use a single scale to decouple the Higgsinos. 



The general approach to solving the RGEs is as 
follows: 

— Begin with weak scale gauge couplings and quark 
masses. Rotate the Yukawas to the current basis 
using four rotation matrices chosen to make sure 
the KM matrix is correct. 

— Run the gauge couplings and Yukawas up to the 
high scale. 

— Input the high scale ansatz which can depend on 
the Yukawa flavour stucture, 

— Run the RGEs down in the presence of thresholds. 
As various thresholds are passed, remove the par- 
ticles from the theory, with special attention given 
to the matching at mn as mentioned above. 

Since the location of the thresholds is not known at the 
beginning of the process, it is necessary to begin with 
an estimate of their location and iterate this procedure 
until the solution converges to the required accuracy. 

In this manner, it is possible to obtain the whole 
set of couplings at a scale appropriate to the prob- 
lem - for example Q = m^^ in the case of stop decay. 
The inputs are: weak scale gauge couplings and squark 
masses; and high scale SUSY parameters governed by 
a particular model. 



3 Additional Couplings when SUSY is 
broken 

When SUSY is broken by the decoupling of one or 
more of the SUSY particles, there are many couplings 
which, although equal in the supersymmetric limit, 
have different RGEs. This causes the couplings them- 
selves to become different. As a simple example we 
can compare the quark-quark-gauge coupling to the 
squark-quark-gaugino coupling. These couplings are 
equal when SUSY is exact, but can be conceptually 
and numerically different below the scale where some 
sparticles have decoupled from the theory. 

To one loop, the running of the gauge coupling does 
not depend on any other couplings and the only effect 
of thresholds is to change the numerical factor in the 
RGE. On the other hand, the running of the gaugino 
coupling contains many extra terms, some depending 
on various gaugino couplings and others depending on 
Yukawa matrices. 

The extra Yukawa terms in the right-handed up- 
type gaugino coupling running are: 

2. -3 [s^, + c'0h] 0Mfu)Jkif:i'')h{9U 

4. -ejfiQXf^)l{g'Q)u{f:''% 

5. 2^^0s,(5:j.fc(/r)r;(/r)y (8) 

In the exact SUSY limit, the extra Yukawa terms can- 
cel and the gaugino couplings are equal to the gauge 
coupling so that the RGEs are identical. When SUSY 
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Fig. 1. The evolution of tlie U(l) gauge coupling to 1-loop 
order (dotted curve) and 2- loop order (dot-dashed curve). 
Also the U(l) gauge coupling with thresholds (solid). The 
upper line (dashed) is the U(l) gaugino coupling to a right- 
handed up-type squark (gi). Vertical dotted lines indicate 
the presence of thresholds. 



is broken this is no longer the case. If the heavy Higgs 
particles decouple first, for example, 9h = and sev- 
eral terms disappear from the running. This means 
that the cancellation of the Yukawa terms will no longer 
be exact, leading to a non-zero overall Yukawa factor 
in the RGE. The coupling can develop non-zero off- 
diagonal flavour changing terms which are not present 
in the gauge coupling. Thus, if the squarks and gaug- 
inos are lighter than mn the RGE for this coupling 
will have a Yukawa component at a scale where it is 
still important for phenomenology. 

In the U(l) gaugino running, the decoupling of 
other gaugino terms can also result in a different value 
for the diagonal elements, as shown in Figure [T] Here 
we show the 1-loop and 2-loop running of gi in the 
dotted curve and dot-dashed curve respectively. The 
solid line shows the running of the gi coupling when 
thresholds are switched on. It is clear that aX Q = mt 
the effect of thresholds on the coupling is comparable 
to the difference between the 1- and 2-loop lines. The 
dashed line running upwards from the gluino thresh- 
old is the g'y^^ coupling which is important until the 
up-right squarks decouple, shown by the central dot- 
ted line. We can see that the gaugino coupling at this 
point is very different from its SM counterpart. 

There are many terms in the full threshold RGEs 
which can provide an unexpected dependence like the 
Yukawa terms in the gaugino coupling RGE. Another 
such example is in the running of the CP-conserving 
gaugino mass, M2 [3, for which 



c.c. 



iA7Tf^3 2M2Sc{~eh 

It is clear that this term does not contribute in the 
SUSY limit (when all 6* = 1) but can appear in the 
RGE if > |^|. 

We have obtained the full system of threshold RGEs 
for the gauge couplings and Yukawas, in addition to 



Table 1. Comparison of the calculated width 
(xlO~^ GeV) for ii cZi using various methods. 

In all cases SSB thresholds were not used and the 
GUT scale inputs were the mSUGRA parameters 
mo = 300 GeV, mi/2 = 250 GeV, ^0 = -1000 GeV, 
tan/3 — 10, sign{jj,) = -1-1, mt = 172 GeV which led to 
mzi = 100 GeV and m^^ ~ 183 GeV. 



Method 



Width 



Hikasa-Kobayashi ~ 18 

1- loop (all thresholds at mn) 1.04 

2- loop (all thresholds at niH) 1.15 
2-loop (realistic thresholds) 1.19 
2-loop (realistic thresholds and tilde terms) 1.31 



threshold RGEs for fj,, gaugino masses, soft masses and 
trilinear SSB parameters [9] and plan to incorporate 
our code into ISAJET (TU]. 



4 Flavour changing top squark decay 

The two body decay of the lighter stop (ii) into a 
charm quark and the lightest neutralino (Zi) occurs 
at the one loop level. This decay becomes important if 
tree level two body decays are kinematically forbidden, 
although it may still compete with three body decays 
TI]. 

Two body flavour changing stop decay was first 
studied by Hikasa and Kobayashi [5], where the off- 
diagonal elements of the up-squark SSB matrix, specif- 
ically i{L,i?} — cl mixing elements, were estimated un- 
der the approximation that the RGEs could be inte- 
grated using a single step and assuming the LSP was 
a photincH. The partial decay width can then be read- 
ily obtained. Although their estimate was adequate for 
the very light stop case they were interested in, it is 
not valid for m^^ in the range of interest today. 

To this end, the RG method was used to calculate 
the mixing numerically. This removes the need for the 
one-step approximation and also allows us to obtain 
the full weak scale couplings for use in the decay cal- 
culation. The decay width was re-derived from the La- 
grangian, keeping the difference between sparticle and 
SM couplings, which allows for the possibility of ex- 
tra flavour changing terms coming from sources other 
than the up-squark mass matrix. 

Table [T] compares this partial width calculated us- 
ing the various methods. The Hikasa-Kobayashi result 
clearly overestimates the width, illustrating the im- 
portance of our improvement. In a random sample of 
10,000 mSUGRA points, over 300 were found with a 
stop mass above 100 GeV and with kinematically for- 
bidden two body decays. In these cases the one-step 
approximation overestimated the width by between a 
factor of 15-25 although in a few cases the difference 
was as large as a factor of 35. 



The Hikasa-Kobayashi result was later modified to al- 
low for arbitrary composition of the LSP. 



Colliders - SUSY Phenomenology 



Contributed Talk 



Looking at the rest of the results in the table, we 
can see that introducing thresholds to the equations 
produces a result which, although less than the dif- 
ference between the 1-loop and 2-loop running, is of 
similar order. This indicates that the thresholds are 
important for claiming true 2-loop accuracy. The in- 
troduction of tilde terms produces a sizeable change, 
which shows that the difference between sparticle and 
SM couplings is also important. 

Clearly the single step integration may give a qual- 
itatively misleading value of the ii decay branching ra- 
tios in the event there are competing modes. The other 
effects are also important for obtaining a quantitative 
projection of stop decay patterns. 



5 Conclusions 

We have constructed a closed system of RG equations 
including threshold decoupling for sparticles and heavy 
Higgs particles. This includes SSB trilinear couplings 
and masses in addition to the gauge and Yukawa cou- 
plings. A code has been developed to solve the RGEs 
which will be incorporated into the Isajet event gen- 
erator. 

When SUSY is broken, the RGEs for sparticle cou- 
plings can become different from the RGEs for their 
SM counterparts, i.e. the gauge and Yukawa couphngs. 
This not only results in different diagonal elements at 
the weak scale, but many sparticle couplings also de- 
velop additional flavour ofl^-diagonal terms. 

Our consideration of the flavour changing decay of 
the stop finds the partial width to have been over- 
estimated by a factor of 15 to 25 as a result of the 
one-step approximation. We also find that threshold 
effects in both the standard couplings and the so-called 
tilde terms are important for more quantitative flavour 
violation calculations. 
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A 4-component Yukawa RGE 

The Lagrangian in ([T]) can be recast into 4-component 
form by combining spinors into 4-component Majo- 
rana fermions given by = (V'-A/i "CV'm) ; ^-nd Dirac 
fermions, 'Poi = (V'l, -CV'^Jr 

^\x]^-WM,PAlk^a + \x]t^M,PLnik^i 

(10) 



where the scalar 4>a is now complex. 

Using this new Lagrangian, the RGEs for t/^" , U^°'^ 
V^W^X^" and X"^"- can be obtained from In the 
case of the MSSM, the Lagrangian does not contain 
any terms where the Higgs fields {K) connect Dirac 
and Majorana fermions, so and vanish. Here 
h can stand for any of the fields h, 7i, or since 
the general form of the RGE is the same for all Higgs 
fields. 

The RGE for [/^'*, which we show as an example, 

is 



(47r) 



\h 



dt 



1 — loop 



-U"Tr {(U"^U1'' + U2''tu2&) 

+ i(X'bt^lh^y^2hf^2. 



W^'^W'')] 
- V^'X^'^tw''] 



2h^ 



+U2''Tr {(U2''tui'' + U^'^tui'') 
+ i(X2^txi'^+X2''tx") 



(U) 



where the final term has been constructed by sepa- 
rating C2{F) into separate terms for left-handed and 
right-handed fields. The other RGEs will appear in [S]. 
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